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We discuss the present collective flow signals for the phase transition to the quark-gluon plasma 
(QGP) and the collective flow as a barometer for the equation of state (EoS). We emphasize the 
importance of the flow excitation function from 1 to 50A GeV: here the hydrodynamic model has 
predicted the collapse of the vi-flow at ^ lOA GeV and of the V2-flow at ~ 40A GeV. In the latter 
case, this has recently been observed by the NA49 collaboration. Since hadronic rescattering 
models predict much larger flow than observed at this energy, we interpret this observation as 
potential evidence for a first order phase transition at high baryon density p^. 
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Figure 1: The phase diagram with the critical end point at Hb « 400 MeV, T « 160 MeV, predicted by 
Lattice QCD calculations. For different bombarding energies, the time evolution in the T — /X^-plane of a 
central cell in UrQMD calculations [nOp is depicted, (from Bratkovskaya et al.) [@]. 



1. The QCD phase diagram 

The phase diagram predicted by lattice QCD calculations [[T], ^ (Fig. [T]) shows a cross over 
for vanishing or small chemical potentials {j.^, but no first-order phase transition to the quark- 
gluon plasma (QGP). This region may be accessible at full RHIC energy. In contrast, at lower 
SPS and RHIC energies {^/s 4 — 12A GeV) and in the fragmentation region of RHIC, y ^3 — 5 
[^, a first-order phase transition is expected with a critical baryochemical potential of [|l], ^ 
400 lb 50 MeV and a critical temperature of w 150 — 160 MeV. This first-order phase 
transition is expected to occur at finite strangeness [^. 

A comparison of the QCD predictions of the thermodynamic parameters T and /ig with the re- 
sults from the UrQMD transport model [|6|, ^ in the central overlap regime of Au-i-Au collisions [§| 
are shown in Figure |l|. The 'experimental' chemical freeze-out parameters - determined from fits 
to the experimental yields - are shown by full dots with errorbars and taken from Ref. [^. The tem- 
perature T and chemical potentials /ig, denoted by triangular and quadratic symbols (time-ordered 
in vertical sequence), are taken from UrQMD transport calculations in central Au-i-Au (Pb-i-Pb) 



collisions at RHIC [10| as a function of the reaction time (separated by 1 fm/c steps from top to 
bottom). Full symbols denote configurations in approximate pressure equilibrium in longitudinal 
and transverse direction, while open symbols denote nonequilibrium configurations and correspond 
to T parameters extracted from the transverse momentum distributions. 

The transport calculations during the nonequilibrium phase (open symbols) show much higher 
temperatures (or energy densities) than the 'experimental' chemical freeze-out configurations at 
all bombarding energies (> llA GeV). These numbers exceed the critical point of (2-1-1) flavor 
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lattice QCD calculations by the Bielefeld-Swansea-collaboration [Q] (large open circle) and by the 
Wuppertal-Budapest-coUaboration [|l]] (open square; the star denotes earlier results from [p). The 
energy density at Hc, Tc is of the order of 1 GeV/fm^. At RHIC energies, when the temperature 
drops during the expansion phase of the 'hot fireball' a cross over is expected at midrapidity. Using 
the statistical model analysis by the BRAHMS collaboration based on measured antibaryon to 
baryon ratios [ 11 ] for different rapidity intervals at RHIC energies, the baryochemical potential pLs 
has been obtained. At midrapidity, one observes pLs — 0, whereas at forward rapidities plb increases 
up to jlB — 130 MeV at J = 3. Thus, only a forward rapidity measurement (j 4 — 5) at RHIC 
will allow to probe large /ig. A unique opportunity to reach higher chemical potentials and the 
first-order phase transition region at midrapidity is offered by the STAR and PHENIX detectors 
at RHIC in the high-/x-RHIC-running at ^ = 4 - 12A GeV. For first results see Ref. [|12|]. The 
International FAIR Facility at GSI will offer a research program fully devoted to this topic in the 
next decade. 



1.1 Flow Effects from Hydrodynamics 



Early in the 70th, hydrodynamic flow and shock formation have been proposed [|T^, |T4|] as 
the key mechanism for the creation of hot and dense matter in relativistic heavy-ion collisions [p^]. 
Though, the full three-dimensional hydrodynamical flow problem is much more complicated than 



the one-dimensional Landau model []16|]. The 3-dimensional compression and expansion dynamics 
yields complex triple differential cross sections which provide quite accurate spectroscopic han- 
dles on the EoS. Differential barometers for the properties of compressed, dense matter from SIS 
to RHIC are the bounce-off, v\{pt) (i.e., the strength of the directed flow in the reaction plane), 
the squeeze-out, V2{pt) (the strength of the second moment of the azimuthal particle emission 



distribution) [ 13, h4L O, h8L h9^, 20, 21], and the antiflow [17, 18, h9^, 20, 21] (third flow com 



ponent [22, ^]). It has been shown [14, 17, 18, [I^, 2^ 21] that the disappearance or so-called 
collapse of flow is a direct result of a first-order phase transition. 



To determine these different barometers, several hydrodynamic models [ |24[ ] have been used 
in the past, starting with the one-fluid ideal hydrodynamic approach. It is known that this model 
predicts far too large flow effects so that viscous fluid models have been developed [p5, 26, 27] to 



obtain a better description of the dynamics. In parallel, so-called three-fluid models, which distin- 



guish between projectile, target and the fireball fluid, have been considered [28]. Here viscosity 
effects do not appear inside the individual fluids, but only between different fluids. One aim is to 



obtain a reliable, three-dimensional, relativistic three-fluid model including viscosity [26, 27]. 

Though flow can be described very elegantly in hydrodynamics, one should consider micro- 
scopic multicomponent (pre-)hadron transport theory, e.g. models like qMD [p9|], IQMD [^], 
UrQMD [g 0], or HSD to control models for viscous hydrodynamics and to gain background 
models to subtract interesting non-hadronic effects from data. If hydrodynamics with and without 
quark matter EoS and hadronic transport models without quark matter - but with strings - are com- 
pared to data, can we learn whether quark matter has been formed? What degree of equilibration 
has been reached? What does the EoS look like? How are the particle properties, self-energies, 
cross sections changed? 
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Figure 2: Sideward flow px of (left) K, A and p's at 6A GeV as measured by E895 in semi-central collisions 



at the AGS and (right) for p and A compared to UrQMD 1.1 calculations for < 7 fm | |39 | 
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Figure 3: Prediction of the directed flow from ideal hydrodynamics with a QGP phase (open symbols) and 
from the Quark Gluon String Model without QGP phase (full symbols) [|2[ . 



1.2 Evidence for a first-order phase transition from AGS and SPS 



The formation and distribution of many hadronic particles at AGS and SPS is quite well de- 



scribed by microscopic (pre-)hadronic transport models [32]. Additionally, flow data are described 



reasonably well up to AGS energies [22, 33, 34, 39, 40, 41 1, if a nuclear potential has been included 
for the low energy regime. 

However, since ideal hydrodynamical calculations predict far too much flow at these ener- 
gies [^], viscosity effects have to be taken into account. While the directed flow Px/m measure- 
ment of the E895 collaboration shows that the p and A data are reproduced reasonably well [39, ?^], 



4 



Collapse of Flow 



Horst Stocker 




Figure 4: The time evolution of directed flow px/N as a function of rapidity for Au+Au collisions at 8A GeV 
in the one-fluid model for (left) a hadronic EoS without phase transition and (right) an EoS including a first- 
order phase transition to the QGP [from Brachmann] . 



ideal hydrodynamical calculations yield factors of two higher values for the sideward flow at 
SIS M\ and AGS. 



However, the appearance of a so-called "third flow component" [ ]22p or "antiflow" [ ]43P in cen- 
tral collisions (cf. Fig. ^ is predicted in ideal hydrodynamics, though only if the matter undergoes 
a first order phase transition to the QGP. It implies that around midrapidity the directed flow, Pxiy), 
of protons develops a negative slope. Such an exotic "antiflow" (negative slope) wiggle in the pro- 
ton flow v\{y) does not appear for a hadronic EoS without QGP phase transition at intermediate 
energies. For high energies see disussion in References [^, Just as the microscopic transport 
theory (Fig. ^r.h.s.) and as the data (Fig. ^l.h.s.), the ideal hydrodynamic time evolution of the 
directed flow, Px/N, for the purely hadronic EoS (Fig. ^ l.h.s.) does show a clean linear increase of 
Px{y)- However, it can be seen that for an EoS including a first order phase transition to the QGP 
(Fig. ^r.h.s.) that the proton flow vi ~ Px/pt collapses around midrapidity. This is explained by 
an antiflow component of protons that develops when the expansion from the plasma sets in p6|]. 

Even negative values of d{px/N)/dy calculated from ideal hydrodynamics (Fig. ||) show up 
between 8 and 20A GeV. An increase up to positive values is predicted with increasing energy. But, 
the hydro calculations suggest this "softest point collapse" is at Eiai, ~ GeV. This predicted 
minimum of the proton flow has not been verified by the AGS data! However, a collapse of the 
directed proton flow at Eiah ~ 30A GeV (Fig. ^ is verified by a linear extrapolation of the AGS 
data. 



5 



Collapse of Flow 



Horst Stocker 




Figure 5: The proton dp^/dy-slope data measured by SIS and AGS compared to a one-fluid hydrodynamical 
calculation. A linear extrapolation of the AGS data indicates a collapse of flow at Ej^,!, w 30A GeV (see also 
Ref. Q). The point at 40A GeV is calculated using the NA49 central data (cf. Alt et al.) [|8]]. 

This prediction has recently been supported by the low energy 40A GeV SPS data of the 



NA49 collaboration (cf. Figs. | and|7|). In contrast to the AGS data as well as to the UrQMD 
calculations involving no phase transition (Figs. ^ and |7|), the first proton "antiflow" around mid- 
rapidity is clearly visible in these data. 

Thus, a first order phase transition to the baryon rich QGP is most likely observed at bombard- 
ing energies of 30 — 40A GeV; e.g. the first order phase transition line in the r-jLt^-diagram has 
been crossed (cf. Fig. |l|). In this energy region, the new FAIR- facility at GSI will operate. It can 
be expected that the baryon flow collapses and other first order QGP phase transition signals can be 
studied soon at the lowest SPS energies as well as at fragmentation region j > 4 — 5 for the RHIC 
and LHC collider energies. At high jig, these experiments will enable a detailed study of the first 
order phase transition as well as of the properties of the baryon rich QGP. 

2. More evidence for a first-order phase transition at highest net baryon densities 

Microscopic transport models, at SIS energies, reproduce the data on the excitation function 
of the proton elliptic flow V2 quite well. The data seem to be described well by a soft, momentum- 



dependent EoS ||47|,^. 

Below ^ 5A GeV, the observed proton flow V2 is smaller than zero, which corresponds to the 
squeeze-out predicted by hydrodynamics long ago [13, 17, 1^, 19, 20, 



From the AGS data, a transition from squeeze-out to in-plane flow in the midrapidity region 
can be seen (Fig. |8|). In accord to the transport caluclations (UrQMD calculations in Fig. |8| [39]; 



for HSD results see [4C, 41]), the proton V2 at 4 — 5A GeV changes its sign. Hadronic transport 



simulations predict a smooth increase of the flow V2 at higher energies (10— 160A GeV). The 
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Figure 6: (Color online)Directed flow of protons (left) and pions (right) in Pb+Pb collisions at fiiab = 
40A GeV with pr <2 GeV/c. UrQMD calculations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method (circles), cumulant method of order 2 (squares) 
and cumulant method of order 3 (triangles) are depicted. The 12.5% most central collisions are labeled 
as central, the centrality 12.5% -33.5% as mid-central and 33.5% -100% as peripheral. For the model 
calculations the corresponding impact parameters of < 3.4 fm for central, = 5 — 9 fm for mid-central and 
^7 = 9 — 15 fm for peripheral collisions have been used (from Petersen et. al. [S9|). 



160A GeV data of the NA49 collaboration indicate that this smooth increase proceeds as predicted 
between AGS and SPS. For midcentral and peripheral protons at 40A GeV (cf. Ref. [|3^, 
UrQMD calculations without phase transition give a considerable 3% V2 flow. 

Contrary, the recent NA49 data at 40A GeV (see Ref. [Q, ^ (cf. Figs. | and [T^) show a 



sudden collapse of the proton flow for midcentral collisions. At 40A GeV this collapse of V2 for 
protons around midrapidity is very pronounced while it is not observed at 16QA GeV 

Another evidence for the hypothesis of the observation of a first-order phase transition to 



QCD is the dramatic collapse of the flow vi also observed by NA49 [38], again around 40A GeV, 
where the collapse of V2 has been observed. This is the highest energy at which a first-order phase 
transition can be reached at central rapidities of relativistic heavy-ion collisions (cf. Ref. [|l], ^ 
and Fig. [IJ). Therefore one may conclude that a first-order phase transition at the highest baryon 
densities accessible in nature has been seen at these energies in Pb-i-Pb collisions. As shown in 



Ref. []50|], the elliptic flow clearly distinguishes between a first-order phase transition and a cross 
over. 



3. Summary 

Evidence for a first-order phase transition in baryon-rich dense matter is recently presented 
by the collapse of both, vi- and V2-collective flow of protons from the Pb-nPb collisions at 40A GeV 
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Figure 7: (Color online)Directed flow of protons (left) and pions (right) in Pb+Pb collisions at fiiab ~ 
160A GeV with pt <2 GeV/c. UrQMD calculations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method (circles), cumulant method of order 2 (squares) 
and cumulant method of order 3 (triangles) are depicted. The 12.5% most central collisions are labeled 
as central, the centrality 12.5% -33.5% as mid-central and 33.5% -100% as peripheral. For the model 
calculations the corresponding impact parameters of < 3.4 fm for central, = 5 — 9 fm for mid-central and 
= 9 — 15 fm for peripheral collisions have been used (from Petersen et. al. [S9|). 



of the NA49 collaboration. It will soon be possible to study the nature of this transition and the 
properties of the QGP at the high-jU/low energy and at the forward fragmentation region at RHIC 
and at the future GSI facihty FAIR. 

This first-order phase transition occurs according to lattice QCD results ^ ^ for chemical 
potentials above 400 MeV. Since the elliptic flow clearly distinguishes between a first-order phase 
transition and a cross over ||50|], the observed collapse of flow, as predicted in Ref. [13, 14], is a 
clear signal for a first-order phase transition at the highest baryon densities. Calculations from ideal 
hydrodynamics including additional fluctuations predict an increase of 50% for fluctuations of 
the flow; however transport models predict an increase by a factor of 2 and 3 [52]. The viscosity 
coefficient of QGP might experimentally be determined from these fluctuations. 

We predict that the collapse of the proton flow analogous to the 40A GeV data will be seen in 
the second-generation experiments at RHIC and FAIR. 
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Figure 8: (Color online) The calculated energy excitation function of elliptic flow of protons in 
Au+Au/Pb+Pb collisions in mid-central collisions (b=5-9 fm) with \y\ < 0.1(full line). This curve is com- 
pared to data from different experiments for mid-central collisions. For E895 [^[^, FOPI and NA49 
there is the elliptic flow of protons. The dotted line in the low energy regime depicts UrQMD calcula- 
tions with included nuclear potential (from Petersen et. al. |^]). 
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